Does the D / D + production asymmetry decrease at large xp? 
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Abstract 

o 

fvi . We have applied the meson cloud model (MCM) to calculate the asymmetries 

^ j in D and D s meson production in high energy £~-nucleus and 7r _ -nucleus 

collisions. We find a good agreement with recent data. Our results suggest 
that the asymmetries may decrease at large xf- 
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Several experiments have reported [|l}|| a significant difference between the xf (Feynman 
momentum) dependence of leading and nonleading charmed mesons (D). Recent data taken 
by the WA89 collaboration [|5| with a £~ beam have not only confirmed the asymmetry and 
the leading particle effect but have also observed this effect in D s and A c production. An 
interesting feature of the WA89 data is that they suggest, inspite of very large error bars and 
poor statistics, that the asymmetry decreases at very large xf- More recently, preliminary 
data from the SELEX collaboration M, in contrast to all previous experiments, seem to 
indicate that the asymmetry is smaller and follows a more complicated pattern. 

Very soon after the appearence of the first data, it became clear that it was not possible to 
understand them only with usual perturbative QCD or with the string fragmentation model 
contained in PYTHIA. Alternative models have been advanced [|7|||. All these models obtain 
a reasonable description of data, but none of them predicts a decrease in the asymmetry. 

The purpose of this letter is to show that in the meson cloud model (MCM) || we can 
reproduce data and accomodate a possible decrease of the asymmetry. The MCM has been 
very successful in the study of hadron structure || and of particle production in high energy 
soft hadron collisions [p~0] , |TT]| . It has been extended to the charm sector ||12|| . 

In the MCM we assume that quantum fluctuations in the projectile play an important 
role. Both the £~ and 7r~ may be decomposed in a series of Fock states. This series has 
also been discussed in Ref. ||, where, for example, fluctuations of the type \tt~) = \udcc) 
and |£~) = \ddscc) were considered. In the MCM we write the Fock decomposition in terms 
of the equivalent hadronic states |7r~) = \D°*D~) and |£~) = |S°D~). This expansion 
contains the "bare" terms (without cloud fluctuations), light states and states containing 
the produced charmed meson (D or D s ). The latter are, of course very much suppressed but 
they will be responsible for asymmetries. The "bare" states occur with a higher probability 
and are responsible for the bulk of charm meson production at low and medium momentum 
(xf < 0.4), including, for example the perturbative QCD contribution. The cloud states 
are less frequent fluctuations and contribute to D production in the ways described below. 
More precisely we shall assume that: 

|E-) = Z[|Eo>+... + |S2£>->+|E°£>7>] (1) 

\n-) = Z' i\n° ) + ... + \D°*D-)] (2) 

where Z and Z' are normalization constants, |Eq) and |7Tg) are the "bare" sigma and pion 
and the "dots" denote possible cloud states \MB) in the E~ and \MM') or \BB ) in the n~ . 
The relative normalization of these states is fixed once the cloud parameters are fixed. We 
shall first study the reactions induced by the cloud component of the E~. This projectile 
baryon is thus regarded as being a sum of virtual meson (M)-baryon (B) pairs and a E~- 
proton reaction can thus be viewed as a reaction between the "constituent" mesons and 
baryons of the E~ with the target proton. 

With a E~ beam the possible reaction mechanisms for D~ meson production at large xf 
and small pt (the soft regime) are illustrated in Fig. 1. 
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FIG. 1. E p collision in which the projectile is in a \MB) state. Figs, la) and lb) show the 
"indirect" D^ (D^ 1 ) production and lc) the "direct" D~ (Dj) production. 



In Fig. la the baryon just "flies through" , whereas the corresponding meson interacts 
inelastically producing a D meson in the final state. In Fig. lb the meson just "flies 
through" , whereas the corresponding baryon interacts inelastically producing a D meson in 
the final state. In Fig. lc the meson in the cloud is already a D~ (or D~) which escapes 
(similar considerations hold for D~ production with a ir~ beam). This last mechanism is 
responsible for generating asymmetries. We shall refer to the first two processes as "indirect 
production" (J) and to the last one as "direct production" (D). The first two are calculated 
with convolution formulas whereas the last one is given basically by the meson momentum 
distribution in the cloud initial \MB) state. Direct production has been widely used in the 
context of the MCM and applied to study n, A ++ and ir° production [|10| . Indirect meson 



production has been considered previously in |TT 



Inside the baryon, in the \MB) state, the meson and baryon have fractional momentum 
yu and ys with distributions called Jmimb^Vm) and /b/mb(|/b) respectively (we shall use for 
them the short notation f M and /#). Of course, by momentum conservation, y M + y B = 1 
and these distributions are related by /m(2/) = /b(1 — y) HHHl- The "splitting function" 



fhi{y) represents the probability density to find a meson with momentum fraction y of the 
total cloud state \MB). With f M and fg we can compute the differential cross section for 
production of D mesons, which, in the reaction S^p — > DX, is given by: 
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where $o and <!>/ refer respectively to "bare" and indirect contributions to D meson produc- 
tion and xf is the fractional longitudinal momentum of the outgoing meson. $£> represents 
the direct process depicted in Fig. lc and is given by |lO|Jlli : 
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f D (x F ) a L 



(4) 



where fu = fD-/s°D- an d u^ is the total pS° cross section. An analogous expression can 
be written for the reaction -K~p — > DX. 



Using (||), we can compute the cross sections and also the leading (D )/nonleading(D + ) 
asymmetry: 



da D {x F ) _ da D (x F ) 



A(x F ) 
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(5) 

where the last line follows from assuming $f = $f = $f\ This last assumption is 
made just for the sake of simplicity. In reality (and also in the calculations performed 
in [13[]) these contributions are not equal and their difference is an additional source of 
asymmetry, which in some cases is not negligible. Since the "bare" states do not give origin 
to D~ / D + asymmetries (they represent mostly perturbative QCD contributions which rarely 
leave quark pairs in the appropriate kinematic region), we have made use of $^ + = $f ~ = 
Qq. The denominator of the above expression can be replaced by a parametrization of the 
experimental data: 

$ T = ob[(l-a*r~ + (l-x F ) n+ ] (6) 

where n~ and n + are powers used by the different collaborations to fit their data and 
a ~ 4 — 7 fib as suggested by the data analysis performed in [0-0]. 
Inserting (Q) and (||) into (||) the asymmetry becomes: 

A(x F ) = ^ — ( ^H -^ (7) 

K ' cto x F [ (i - x F y- + (i - x F ) n+ ] v; 

The behavior of (|7]) is controlled by fn{x F )- In the recent works with the MCM one 
finds two forms for the splitting functions. One comes from the evaluation of the relevant 
Feynman diagrams (Sullivan process) [f|,li] and the other comes from a light cone ansatz for 
the cloud state wave function |14[] . We shall compute the asymmetries with both of them. 
The light cone splitting function is given by 



fM{v) = T, / iidsliMv,ki)l s 



l ;HWy{l-y)eM-^k) (8) 



4vr 2 " v " y " v 4a 2 

where ij)\\>(y,\<L\) is the light cone wave function of the Fock state containing a meson M 
(baryon B), with longitudinal momentum fraction y (or 1 — y), transverse momentum k^ 

(— kj_) and helicity A (A'); M.\i — — "*■ H — \_ B is the invariant mass of the meson (M) 

- baryon (B) system for large longitudinal momenta, rriyi and mg are their masses, if is a 
normalization constant and a is the width of the heavy meson-baryon state. 



Substituting @ into (0) and remembering that in our case y = x F , we can write our 
final expression for the asymmetry: 

A 1 (x F )=N 1 {1 ~ XF) 



[{i-x F y- + (\-x F y + ] 

where Ni = H 2 a 2 <r s /4 tx a . 

A striking feature of (||) is that in the limit x F — *■ 1 it goes to zero regardless of the choice 
of two parameters N\ and a. This happens because when the leading meson (responsible for 
the asymmetry) has momentum one its accompanying cloud baryon must have momentum 
zero, being strongly virtual, thereby increasing the invariant mass of the cloud state and 
forcing this cloud configuration to have zero probability. For the pion beam, we obtain 
an analogous expression with the replacements m= — > m D o*, o~ — > o D . In this case 
the parameters H and a may assume different values. Apart from numerical changes, the 
qualitative behavior of A\ remains the same, also for large x F . 

We now write the splitting function in the Sullivan process approach. The fractional 
momentum distribution of a pseudoscalar meson M in the state \MB' > is given by p.|TT 



9mbb> f tmax j. [—t + ( m B> -m B ) 2 } 



167T 2 J-oc \t — mi- 1 - 



x F 2 MBB ,{t) (10) 

where t and tum are the four momentum square and the mass of the meson in the cloud state 
and t m ax = Tn% y — m 2 B , y/(l—y) is the maximum t, with tub and tub 1 respectively the B and 
B' masses. Following a phenomenological approach, we use for the baryon-meson-baryon 
form factor Fmbb 1 , the exponential form: 

F MBB >(t) = exp fc^l (11) 

V 1Y MBB' J 

where Ambb' is the form factor cut-off parameter. 

Considering the particular case where B = £~, B' = 2° and M = D~, we insert ( |ToD 
into (0) to obtain the final expression for the asymmetry in our second approach: 

A 2 (x F ) 



[(l-x F ) n - + (l-x F ) n+ 

tmax [— t + (jJJn — 7B£) 2 1 



x L dt - V-W' ' Flss{t) (12) 

where N 2 = ^ss cr"/16 7r a - 

For the pion beam, we need also the Sullivan splitting function of the state |7r~ > — > 
\D°*D~). In this state, the D meson momentum distribution (which was computed in Ref. 
fl3||) turns out to be identical to (pTOj) except for the bracket in the numerator which takes 



the form |T5fl , [— t + ((m 2 — m D a* — t)/2m D o*) ], and for trivial changes in the definitions, 



i- e -> 9mbb> -> 91dd°*i Fmbb' -» F wDD o, and A M bb> -» Kdd°*- Realizing that y = x F in 
the above equations, we can see that in the limit xf — > 1, t ma:r — > — oo and the integral in 
(|I2|) goes to zero. In fact, it vanishes faster than the denominator and therefore A 2 — > 0. 
This behavior does not depend on the cut-off parameter but it depends on the choice of 
the form factor. For a monopole form factor we may obtain asymmetries which grow even 
at very large Xf- Since t controls the off-shellness of the emitted meson, which, in turn, 
is related to the virtuality of the \MB' > state (or \MM' > state in the case of the pion 
beam), the vanishing of A 2 happens for the same physical reason of the vanishing of A x . 

The Dj/Df production asymmetry (with the S~ beam) can be calculated following the 
steps mentioned above and replacing the \E? C D~) state by |E°Z?~). Of course, this implies 
different values for a and N\ (in the light cone approach) and for A and N 2 (in the Sullivan 
approach), but the qualitative discussion (and conclusions) presented above for A\ and A 2 
remain valid. 

Before presenting our numerical results, we emphasize that i) our calculation is based 
on quite general and well established ideas, namely that hadron projectiles fluctuate into 
hadron-hadron (cloud) states and that these states interact with the target; ii) our results 
only depend on two parameters: A\ depends on a and N% and A 2 depends on A and 
A^ 2 . Whereas a and A affect the width and position of the maximum of the momentum 
distribution of the leading meson in the cloud (and consequently of the asymmetry), N x and 
A^ 2 are multiplicative constants which determine the strength of the asymmetry. 

We show in Figs. 2a and 2b, respectively, our results for the asymmetries D~ / D + and 
D~ /Df (for a S~ beam). The results for a ir~ beam are shown in Fig. 3. In all figures solid 
(dashed) lines represent A 2 (Ax). In both approaches we have two parameters which may 
be different from reaction to reaction. They are given in Tables 1 and 2 below: 



Eq. (|) 


a (GeV) 


n + 


n 


N x 


a D /a T {%) 


£-A^ D~X 


0.77 


5.0 


4.5 


400.00 


13.0 


YrA -> D~X 


0.47 


7.0 


4.0 


0.72 x 10 9 


9.0 


ii- A -► D~X 


1.20 


5.0 


3.5 


1.10 


8.0 



Table 1. Parameters used in the asymmetry A\ 
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A (GeV) 


n + 


n 


A^2 


* n /* T {%) 


TrA -► D'X 


2.64 


5.0 


4.5 


2.40 


52.0 


TrA -> D~X 


2.52 


7.0 


4.0 


5.20 


67.0 


7i~ A ^ D~X 


2.88 


5.0 


3.5 


0.82 


20.0 



Table 2. Parameters used in the asymmetry A 2 
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FIG. 2. (2a) Comparison of the MCM asymmetry, Eq. ([/]), with experimental data 
D-/D+; (2b) the same as (2a) for D~/D+. 
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FIG. 3. Comparison of the MCM asymmetry, Eq. (|7]), with experimental data [|i]-[i for D /D + . 

As it can be seen in Figs. 2 and 3 the agreement between the MCM and data is very good. 
In our picture it is simple to understand why the D~ / D + asymmetry peaks at Xp ~ 0.55 
for the £~ beam (Fig. 2a) whereas it peaks at a much larger value xp ~ 0.8 for the ir~ 
beam (Fig. 3). The D~ meson in the S~ beam originates from the |S°Z)~) state and in 
the 7T~ beam it comes from the \D°*D~) state. Since in the meson-meson state the masses 
are closer than they are in the baryon-meson state, the D~ is "faster" inside the pion than 



inside the E~. In ([F]), /d(^f) will peak sooner and die faster for the E~. We emphasize 
that what shifts the peaks of the asymmetries are the masses involved rather than the a (or 
A) parameter. This makes the overall behavior of the asymmetries weakly dependent on 
parameter choices. 

The last column in both tables shows the ratio between the direct and total D~ cross 
sections (where a T = o"o / d%F (1 — xf) h )• We find this quantity quite model dependent. 
This is also a consistency check. In order to treat the cloud as a perturbation we expect this 
ratio to be at the level of 10 — 20%. In two cases, we observe a significant deviation from this 
expectation. This means that, in these cases, we need a large normalization constant for the 
cloud state responsible for direct production (which implies large N 2 ), in order to reproduce 
the observed asymmetry. This is a consequence of neglecting the asymmetry generated by 
indirect production, i.e., in these cases, the approximation <&f = <&f is not a good one. 



A more complete discussion, taking this fact into account is presented in [13|. Here, for the 
sake of the argument, we prefer to keep the calculation simple and keep our fits as they are, 
being implicit that the large A^'s mimic the inclusion of indirect production. 

In conclusion, we have shown that the MCM provides a good understanding of the 
charm production asymmetries in terms of a simple physical picture with few parameters. 
It connects the behavior of the asymmetries at large xp with the charm meson momentum 
distribution within the cloud state. We can explain why we observe asymmetries, why they 
are different for different beams and we are led to the conclusion that they may vanish at 
very large xf- This approach can easily be extendend to charm baryon production and also 
to B production fHf] , 
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